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CONVECTIVE  STABILIZATION  OF  IONOSPHERIC  PLASMA  CLOUDS 


I.  INTRODUCTION 

The  evolution  of  artificial  pla^a  clouds  (e.g. ,  barium)  in  the 
earth’s  ionosphere  continues  to  be  of  interest  to  space  plasma  physicists 
after  more  than  two  decades  of  research.  The  plasma  dynamics  associated 
with  ionospheric  clouds  are  of  interest  since  they  provide  a  diagnostic  of 
the  ambient  ionospheric  environment,  and  also  provide  experimental  data  for 
the  study  of  plasma  instabilities.  Of  particular  interest  is  the  onset  and 
evolution  of  the  E  *  B  gradient  drift  instability  (Simon,  1963;  Hoh,  1963) 
which  is  believed  to  cause  the  gross  structuring  of  plasma  clouds,  i.e., 
field-aligned  strlatlons  (Llnson  and  Workman,  1970).  This  instability  is 
also  believed  to  be  responsible  for,  at  times,  the  structuring  of  the 
ambient,  high  latitude  F  region  lonosj^ere  (Kesklnen  and  Ossakow,  1983). 
The  E  *  B  gradient  drift  instability  is  an  interchange  mode  and  is  driven 
by  a  neutral  wind  or  dc  electric  field  in  an  inhomogeneous,  weakly 
colllsional  plasma.  A  substantial,  amount  of  theoretical  and  computational 
research  has  been  carried  out  to  understand  the  E  *  B  gradient  drift 
instability  and  its  relevance  to  ionospheric  structure  (Volk  and  Haerendel, 
1971;  Perkins  et  al.,  1973:  Zabusky  et  al.,  1973;  Shiau  and  Simon,  1972; 
Perkins  and  Doles,  1975;  Scannapieco  et  al.,  1976;  Qiaturvedi  and  Ossakow, 
1979;  Keskinen  et  al.,  1980;  McDonald  et  al.,  1980,  1981;  Huba  et  al., 
1983;  Sperling,  1983,  198^1;  Overman  et  al.,  1983;  Sperling  and  Classman, 
1985;  Sperling  et  al.,  1984;  Drake  et  al.,  1985). 

The  bulk  of  theoretical  analyses  to  date  have  been  based  on  the  local 

approximation  in  slab  geometry.  The  local  approximation  la  appropriate  for 

unstable  modes  which  have  wavelengths  much  shorter  than  the  scale  length  of 
Manuacript  approved  November  19,  1986. 


the  density  gradient  associated  with  the  cloud  boundary  (i.e.,  kL  >>  1 
where  k  is  the  wavenumber  and  L  is  the  density  gradient  scale  length)  (Huba 
9t  al.,  1983,  Huba  and  Zalesak,  1983).  However,  the  observed  gross 
structuring  of  plasma  clouds  seems  to  suggest  that  the  dominant  modes  have 
kL  <  1  so  that  the  local  approximation  may  not  be  valid.  The  slab  or  one-> 
dimensional  [n  -  ng(x)]  models  of  plasma  clouds  do  not  include  a  number  of 
important  physical  affects  which  appear  in  more  realistic  two-dimensional 
models.  The  polarization  of  large  plasma  clouds  in  20  greatly  reduces  the 
relative  slip  velocity  of  the  cloud  and  the  ambient  neutral  wind  thereby 
weakening  the  E  >«  B  gradient  drift  instability  (Overman  et  al.,  1983; 
Zalesak  and  Kuba,  1984).  Furthermore,  convection  and/or  propagation  of  the 
perturbations  from  the  unstable  backside  to  the  stable  frontside  of  the 
cloud  may  influence  the  overall  stability  of  20  models. 

The  most  detailed  linear  stability  analysis  of  20  plasma  clouds  in  the 
long  wavelength  regime  has  been  based  on  the  waterbag  model  (Overman  et 
al.,  1983).  The  purpose  of  this  paper  is  to  extend  the  analysis  of  Overman 
et  al.  (1983)  by  including  parallel  dynamics,  i.e.,  density  and  potential 
fluctuations  along  the  ambient  magnetic  field.  Recently  it  has  been  shown 
that  the  parallel  effects  can  strongly  Influence  the  linear  stability  of 
ionospheric  plasma  clouds  by  stabilizing  the  short  wavelength  modes 
(Sperling,  1983;  Sperling  et  al.,  1984;  Sperling  and  Classman,  1985;  Drake 
et  al.,  1985).  However,  these  investigations  have  been  based  on  a  1 0  cloud 
model  and  the  need  to  consider  a  20  cloud  model  is  apparent.  Thus,  in  this 
paper  we  study  the  influence  of  parallel  dynamics  on  the  stability  of  20 
ionospheric  plasma  clouds.  In  particular,  as  in  Overman  et  al.  (1983)  we 
consider  a  simple  equilibrium  consisting  of  a  20  cylindrical  waterbag,  but 
in  contrast  to  them  we  consider  three  dimensional  perturbations.  We  find 


that  the  finite  parallel  dynamics  can  dramatically  alter  the  E  *  B  gradient 

drift  instability.  Unlike  the  limit  -  0,  exponentially  growing  global 

eigenmodes  can  exist  when  k^  •  0.  As  k^  is  increased  above  a  threshold 

value,  the  unstable  modes  localize  at  a  finite  angle  away  from  the 

backside;  at  a  point  where  the  diamagnetic  propagation  velocity  (V^j) 

balances  the  convective  flow  velocity  of  the  background  plasma  around  the 

cloud  (Vjj).  We  find  that  the  E  *  B  gradient  drift  instability  is  stable 

when  >  Vjj  so  that  the  cloud  is  no  longer  susceptible  to  large-scale 

structuring.  We  apply  these  results  to  ionospheric  barium  clouds  and 

estimate  that  they  will  cease  structuring  when  L  <  (cT/eB)(M+2)/2V  where 

1  " 

L  is  the  transverse  size  of  the  cloud,  T  ■  T_  +  T.  is  the  total 

X  e  1 

temperature,  M  ■  n^/njj,  and  n^  la  the  cloud  density,  njj  is  the  background 

density,  and  is  the  neutral  wind  velocity.  For  mid-latitude  barium 

releases  at  -  180  km  we  estimate  L  -  160  480  m  which  is  consistent  with 

X 

observations. 

The  organization  of  the  paper  is  as  follows.  In  the  next  section  we 
present  the  assumptions  and  general  equations  used  in  the  analysis.  In 

Section  III  we  derive  the  dispersion  equation  and  in  Section  IV  we  present 

analytical  and  numerical  results.  In  Section  V  we  investigate  the 
possibility  of  exponentially  growing  global  eigenmodes.  Finally,  in 
Section  VI  we  summarize  our  findings  and  discuss  the  application  of  our 
theory  to  the  evolution  of  ionospheric  barium  clouds. 


II.  GENERAL  EQUATIONS  AND  EQUILIBRIUM 


The  general  three  dimensional  equations  for  a  warm  plasma  cloud  In  a 

A  A 

uniform  magnetic  field  B  -  B  e^  and  a  background  neutral  wind  e^ 

(see  Fig.  1)  are  given  by  (Drake  et  al.,  1985) 
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where  n  -  m  v  /ne  is  the  parallel  resistivity,  v  •  +  v  D  .  - 

eee  eeienj_i 

[vin/ilijcTi/eB  is  the  perpendicular  ion  diffusion  coefficient,  is  the 

electronj-ion  collision  frequency  and  Q  and  v  are  the  cyclotron  and 

a  otfi 

neutral  collision  frequencies  of  the  species  a.  Equation  (1)  is  the 

electron  continuity  equation  and  (2)  arises  from  charge  neutrality 
(7*J  ■  0).  We  have  considered  the  electrostatic  limit,  have  assumed 
V  /Q  ,  V.  /fl,  <<  1,  and  have  neglected  the  ion  parallel  diffusion  [d  )  and 
perpendicular  electron  diffusion  [d  ].  We  therefore  assume  that 


3/3t  »  D  ,  3^/3z^, 

ii 


It  is  convenient  to  change  variables  by  defining  a  new  potential 
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Equations  (1)  and  (2)  then  become 
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where  T  -  Tg  +  Tj^. 

We  consider  a  simple  equilibrium  consisting  of  a  cylindrical  waterbag 
of  radius  r^,  and  density  n^  in  a  uniform  background  n^^: 


no(r) 


n  H(r  -  r} 
c  '■  c  ■' 


+  n. 


(6) 


where  H  is  the  Heaviside  function.  The  solution  of  the  equilibrium 
equations  for  this  configuration  are  well  known.  The  drag  between  the 
cloud  and  neutral  wind  polarizes  the  cloud.  The  resulting  potential  is 
given  by  (Sraythe,  1950), 


''n  mT2 


r  <  r 
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where  M  »  n^/n^.  The  potential  causes  the  cloud  to  drift  with  a  uniform 
velocity 


V  *  V  A 

-c  n  M+2  'x' 


(7) 


The  linear  stability  analysis  which  follows  is  most  easily  carried  out 


in  the  frame  of  reference  of  the  moving  cloud.  In  this  frame  the 
potential  is  given  by 


V 
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III.  LINEARIZED  EQUATIONS  AND  DISPERSION  EQUATION. 

We  investigate  the  stability  of  this  equilibrium  by  considering  small 
perturbations  n(r,0,z,t)  and  ♦(r,9,z,t)  around  nQ  and  Since  the 

equilibrium  is  independent  of  z,  we  can  expand  the  perturbations  in  plane 
waves  in  the  z  direction  (n,  ♦  -  exp  ik^^)  without  loss  of  generality.  In 
the  perpendicular  plane  such  a  simple  expansion  is  not  generally  possible 
since  the  equilibrium  depends  on  both  r  and  9.  Equations  (4)  and  (5)  yield 
the  linearized  equations, 
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Equations  (9)  are  solved  separately  in  the  region  r  >  r^,  and  r  <  r^,,  and 

the  solutions  are  then  matched  using  boundary  conditions  which  are  obtained 

from  the  equations  for  r  ■  r^,.  In  the  limit  -  0,  this  procedure  is 

2  - 

straightforward  since  n  -  0  and  7  $  -  0  for  r  •  r^.  In  the  case  k^  *  0, 
the  equations  for  n  and  J  in  the  region  r  •  r^  are  much  more  complicated  so 


further  simplifications  must  be  made.  We  limit  the  calculas-tion  to 
perturbations  for  which  3/39  »  1.  In  this  limit  the  0  dependence 
of  n  and  $  can  be  represented  by  the  eikonal  n  -  exp[iS(0)]  where 

7S  -  iKg(0)  Sg.  (10) 

In  this  limit  the  perturbations  are  strongly  localized  around  the  boundary 
r.  and  decay  exponentially  away  from  this  boundary.  Thus,  for  r  >•  r^  we 
write 


-  (()^  exp[iS(9)]exp[^  k^^[r  •-  r^)] 


(11) 


where  +  and  refer  to  the  region  r  >  r^  and  r  <  r^,  respectively.  The 

specific  range  of  parameters  for  which  the  form  of  i  given  in  (11)  is  valid 

will  be  presented  later.  Similar  expressions  can  be  written  for  n+(r,  0). 
Finally,  we  assume  n  and  5  grow  exponentially  in  time  with  a  growth  rate 
Y.  For  the  present,  Y  must  be  considered  a  local  growth  rate  Y(9). 
Eventually,  we  will  investigate  under  what  circumstances  exponentially 
growing  normal  modes  of  the  cloud  can  exist. 

With  the  form  of  i  (and  n)  given  in  (11),  (9a)  and  (9b)  reduce  to 

algebraic  expressions  for  in  the  regions  r  •  r^,. 
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a  -  £J  Q  /v  V.  >>  1,  D  -  T/m  \>  is  the  parallel  electron  diffusion 
e  i  e  in  .  |e  e  e 

coefficient  and 


Y*  n^  -  o^./eng^. 


The  matching  conditions  for  n^,  and  at  r  -  r^  can  be  derived 
from  (9).  The  radial  E  x  B  convection  of  the  cloud  causes  the  density  n  to 
be  singular  at  the  cloud  boundary  (third  term  in  the  continuity  equation). 
It  is  convenient  therefore  to  separate  out  the  singular  behavior  of  n  so 
that  in  the  region  r  -  r- 


A 

n(r,  0)  -  [ii  5(rrtr  }  +  n^]  exp[iS(9)]. 


The  potential  «  remains  finite  at  the  boundary  so  (9a)  can  be  written  as 


[y  +  k^  +  (c/B)  ikg  (c/B)  ik^n^  ♦  -  0. 


for  r  -  r^,.  Equation  (15)  now  contains  no  singularities  at  r  -  r^. 

However,  is  discontinuous  across  r  -  r^  since  it  is  zero  for  r  <  r^  and 
uc 

finite  for  r  >  r^.  The  parallel  diffusion  coefficient  is  also 

generally  not  continuous  across  the  boundau'y.  These  discontinuities  must 
be  balanced  by  a  corresponding  jump  in  the  potential  ♦  at  the  boundary. 


Thus,  from  (8)  and  (15),  we  find 
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This  condition  can  also  be  derived  from  the  requirement  that  the  tangential 
electric  field  across  the  boundary  of  the  plasma  cloud  be  continuous.  The 
continuity  equation  then  reduces  to 


Y  N  +  (c/B)ik^n  ♦  -  0. 

0  C 


(17) 


The  matching  condition  for  can  be  derived  by  integrating  (9b)  across 
the  boundary. 
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Equations  (13)  and  (16)-(18)  can  then  be  combined  into  a  single  relation 

between  k  and  k  , 
r+  r- 
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Equations  (12)  and  (19)  constitute  a  local  dispersion  equation  for  the 
E  X  B  instability  with  k_  •  0  based  on  a  2D  waterbag  equilibrium. 


IV.  EVALUATION  OF  THE  LOCAL  DISPERSION  EQUATION 

We  first  evaluate  the  local  growth  rate  Y(0)  by  solving  the  dispersion 
equation  given  by  (12)  and  (19).  We  consider  two  limits:  the  cold  plasma 
limit  (T  ■  0)  and  the  warm  plasma  limit  (T  •  0). 

A.  Cold  Plasma  Limit :  T  «  0 

In  the  cold  plasma  limit  (19)  simplifies  to 

2V  V  2 

k  *  (1+M)k  Y  -  [Mkf  COS0  +  (l-^M}(a^  -■  a  )  ^  k  sine]  (20a) 
r+  +  r-i  M+2  '■0  v  /v  +  z 

where 


Y^  -  Y  +  2ik.MV  3in0/(M+2) 
0  n 


(20b) 


Y  -  Y 


(20c) 


since  D  <*  T  ■  0.  Equations  (12c)  and  (12b)  for  k„.  remain  unchanged.  In 

I® 

the  limit  of  small  k^t  (12a)  and  (12b)  yield 


n 


and  the  growth  rate  Y  can  then  be  calculated  from  (20), 


(21a) 


Y«  mY-fcos0  -  i  3in0j 


(21b) 


where  m  -  k.r  is  the  poloidal  mode  number.  The  growth  rate  peaks  at  0  -  0 
0  c 

(the  backside  of  the  plasma  cloud)  and  is  linearly  proportional  to  m.  The 

growth  rate  in  (21)  is  consistent  with  previous  rigorous  and  heuristic 

investigations  of  the  stability  of  circular  waterbag  models  of  plasma 

clouds  (Overman  et  al.,  1983;  Zalesak  and  Huba,  1984).  For  9  •  0,  the 

second  term  on  the  right  side  of  (21b)  causes  the  mode  to  propagate  at  a 

finite  frequency.  This  propagation  results  as  the  fluid  outside  of  the 

cloud  convects  around  the  circular  boundary  and  carries  the  perturbation. 

The  point  9-0  corresponds  to  a  stagnation  point  of  the  flow  so  there  is 

no  finite  frequency  there.  Finally,  from  (21a)  the  assumption  that  the 

modes  are  strongly  localized  around  the  cloud  boundary  requires  m  >>  1  . 

To  obtain  analytic  expressions  for  the  growth  rate  of  modes  with 

k  *  0,  we  consider  only  the  case  9-0  which  corresponds  to  the  most 
z 

unstable  mode.  Later  we  will  present  numerical  solutions  of  the  more 

general  dispersion  relation  in  which  this  restriction  is  relaxed.  For 

simplicity,  we  also  assume  that  v  *  v  >>  v  ,  so  that  a  “  o..  ■  ol 

eenei 

For  9  »  0,  we  have  k  ^  -  k  ■  k  and  -  T  •  Y  and  the  dispersion 

P  +  P"  P  T 

equations  are  given  by, 


Yk  r  -  m  Y-, 
r  c  0’ 


(22a) 


2  2  2  ,2r,  ,  M+2  ''in- 


Vc  “  ”  "  ^  ^  “"y-— J 

1 


(22b) 


with  the  solution 


Y  -  mY  ([m  [m^  +  r) 


2i-I 


i.'-  TH  UTt*  L^!.^  L’Wl'lgSm'WX'^  LI  i^  W  L'.  Wiim'.  'J 


where 


,  1/2, 

^  ‘  “b  ^''c* 


(23b) 


As  is  Increased  from  zero,  over  the  interval  A  <  m  the  growth  rate  is 
roughly  given  by  its  k  •  0  limit.  For 


2  2  2  2 
1  <<  X  /m  «  Q./v,^ 
i  in 


(2iia) 


the  growth  rate  decreases  with  k_  as 


y  -  m  Yq/X. 


(24b) 


1/2 


Over  this  range  k^  *0,1 k  »  k„  so  that  k_  causes  the  mode  to  become 
r  b  z  9  * 

more  localized  near  the  cloud  boundary.  For 


r?.'.vA 

^■WLVV 


,22^ 

b  z  c 


4  M 


2  n: 


(M.2)2 


(25a) 


the  mode  is  stable.  Note  that  (25a)  implies  that  the  lowest  order  poloidal 
modes  are  stabilized  first  as  k_  is  increased.  The  stability  condition 

4t 


given  in  (25b)  is  only  valid  for  << 


In  the  opposite 


limit  >>  V  a 
el  en 


a^/(1+M)  and  the  stability  condition  can  be 
b 


similarly  derived, 


,22. 

b  z  c 


4M 


a 


[l  ♦  (M+l 


i  2 

[D  • 


(25b) 


*  A 
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To  confirm  the  analytic  results  for  T  -  0  and  to  generalize  the 

results  for  0  •  0,  we  present  Figs.  2  and  3  which  are  obtained  by  solving 

(12)  and  (19)  numerically.  In  each  figure  we  plot  (a)  the  growth  rate 
1  /2 

Y/Yq  vs.  X  *  r'eal  frequency  vs.  X  for  the 

parameters  m»1,  v,/n,  -  0.025,  M  -  2,  and  0  -  0  (i.e.,  T  -  0),  and 

several  values  of  0:  (A)  0  -  0®,  (B)  0  »  22.5®,  (C)  0  -  'IS®,  (D)  9  -  67.5®, 

and  (E)  0  -  90®.  In  Fig.  2  we  have  taken  the  limit  v  *  v  >>  v  .  while 

e  en  ei 

in  Fig.  3  we  consider  the  opposite  limit,  v  »  v  .  >>  v 

0  0L  0n 

We  note  the  following  from  Fig.  2  (the  limit  >>  v^j^).  First,  the 
growth  rate  Y  is  a  maximum  when  -  0,  and  decreases  as  9  increases. 

For  9  »  90®  there  is  no  growth.  On  the  other  hand,  the  real  frequency  u 

r 

increases  in  magnitude  as  a  function  of  9  for  X  -  0,  i.e.,  -  0.  These 

points  are  consistent  with  (21).  Second,  for  all  values  of  9,  Y  decreases 

as  X  increases  which  is  consistent  with  (24b),  while  remains  roughly 

constant.  Finally,  the  modes  become  purely  propagating  for  sufficiently 

large  X,  i.e.,  Y  •  0  but  •  0.  The  value  of  X  denoted  by  the  arrow  in 

Fig.  2a  corresponds  to  the  stabilization  point  predicted  by  (25a)  which  is 

in  excellent  agreement  with  the  numerical  results.  We  also  note  that  the 

stabilization  point  is  independent  of  9. 

We  now  discuss  Fig.  3  (the  limit  v  .  >>  v  ).  First,  for  X  -  0  we 

ei  en 

note  that  Y  decreases  as  8  increases,  while  u  increases  in  magnitude:  this 

r 

is  similar  to  the  results  in  Fig.  2.  Second,  in  contrast  to  Fig.  2,  we 

find  that  a  secondary  maxima  occurs  in  Y  at  X  *  36,  and  that  for  finite  k^, 

modes  at  9  •  90®  are  unstable.  This  enhanced  instability  at  finite  values 

of  i3  due  to  the  term  proportional  to  (a^  -  a_)sin  0  in  (20a).  For 

V  >>  V  .  we  have  a  -  a  so  that  this  term  does  not  contribute  to  the 
en  ei 

growth  of  the  mode.  However,  for  v  ,  >>  v  we  note  that  #  a  so  that 

ei  en  - 


14 


this  term  la  finite  and  positive  for  0  <  9  <  ir.  Finally,  the  modea  are 
stable  for  sufficiently  large  X.  The  arrow  in  Fig,  3a  denotes  the 
stability  condition  given  in  (25b)  for  9-0®;  again  this  value  is  in 
excellent  agreement  with  the  numerical  results. 

B.  Warm  Plasma  Limit:  T  •  0 

We  only  analytically  investigate  modes  at  9  -  0  and  again  consider  the 

limit  V  >>  V  , .  Equations  (12)  and  (19)  reduce  to 
en  el 


k  r  Y  -  m(mY-  i  ~  —  k^D  1 

r  c  -  ^0  M+2  2  le-' 


(26a) 


2  2  2  2 
k  r  -  m  X  Y/Yrt 
r  0 


(26b) 


where  k^.^,  -  k^,  Y^  -  Y_  •  Y  +  have  neglected  the  terms 

proportional  to  in  (26b)  for  k^.  This  approximation  is  justified 

for  sufficiently  large  T  (to  be  proven  later). 

For  kg  ■»  0,  the  finite  temperature  corrections  drop  out  of  (26)  and 
the  growth  rate  reduces  to  Y  -  m  Y^  as  in  the  cold  plasma  limit.  As  kg 

increases  the  second  term  on  the  right  side  of  (26a)  becomes  comparable  to 

2  2 
the  driving  term  Y-  when  k,D  -  Y„  v,  /a,  <<  Y-.  At  this  point  k  D  can 
0  z  |e  0  in  i  0  2  |e 

still  be  neglected  compared  with  Y  -  Y^  in  the  definition  of  Y_.  Thus,  the 
eigenvalue  Y  is  given  by 


where 


“'0  ' 

-2  ^  *2)1/2 


(27a) 


ra^  >  X^ 


(27b) 


^  ***  ***  *•*  _• ' ~ 


satisfied,  finite  temperature  effects  stabilize  the  |  >«  B  instability,  and 
the  stability  point  is  given  in  (28b).  In  the  opposite  liiul*" 


•'e"!  '•u.'A 


(29) 


the  cold  plasma  effects  stabilize  the  mode  and  the  stability  point  is  given 

in  (25).  An  expression  similar  to  that  given  in  (28b)  can  again  be  derived 

in  the  limit  u  ,  >>  v  . 

ei  en 

To  verify  our  analytic  results  for  T  •  0,  and  to  extend  these  results 

to  9  •  0,  as  well  as  to  both  v  >>  v  .  and  v  , >>  v  ,  we  present  numerical 

en  ei  ei  en 

solutions  to  the  complete  dispersion  equation  [(12)  and  (19)].  In  Figs.  4 
and  5  we  plot  (a)  the  growth  rate  Y  vs.  X,  and  (b)  the  real  frequency 

w  2 

vs.  X  for  the  parameters  m»1,  v. /fl.  -  0.025,  M  -  2,  /a.  r  Y-  »  5.0  x 

li  X  1 6  D  C  U 

10'^,  and  several  values  of  0;  (A)  9  -  0®,  (B)  9  ■  22.5®,  (C)  9  *  45®, 

(D)  9  -  67.5®,  and  (E)  0  -  90®.  In  Fig.  we  4  consider  the  limit  * 

V  >>  V  , ,  while  in  Fig.  5  we  take  v  *  v  ,  >>  v  .  Thus,  Fig.  4  corres- 
en  ei  e  el  en 

pouds  to  the  warm  temperature  limit  of  Fig.  2  and  Fig.  5  corresponds  to 
that  of  Fig.  3. 


In  Fig.  4  we  note  several  similarities  to  Fig.  2.  The  growth  rate  Y 


decreases  as  9  or  k_  increases  and  the  mode  is  stable  (Y  <  0)  for 

2 


sufficiently  large  k_.  On  the  other  hand,  for  T  •  0  (or  D  *0)  the  mode 

z  ie 


stabilizes  at  a  much  smaller  value  of  k^  (other  parameters  being  equal)  and 


the  stability  point  is  sensitive  to  9.  The  arrow  denotes  the  value  of  X 

for  marginal  stability  (Y  -  0)  based  on  (28b)  for  9  -  0®.  This  value 

(X  -  15.8)  la  in  very  good  agreement  with  the  numerical  results  (X  •  15.5). 

In  Fig.  5  there  are  also  similarities  to  Fig.  3.  For  k^  •  0,  the 
growth  rate  Y  decreases  and  the  real  frequency  ui  increases  in  magnitude 


j  ■.  r: 


E^- 


as  9  increases.  A  secondary  maxima  (or  plateau)  in  Y  occurs  for  a  finite 
value  of  k^.  However,  in  contrast  to  Fig.  2,  the  value  of  k^  at  marginal 
stability  increases  with  9  so  that  the  most  difficult  mode  to  stabilize  is 
at  9  »  90°.  Also,  the  value  of  at  marginal  stability  is  smaller  than 
when  T  »  0. 

It  is  also  illuminating  to  calculate  the  damping  rate  of  the 
instability  in  the  limit  of  very  large  k^.  There  are  two  modes  of  the 
system.  The  first  damps  at  the  ion  diffusion  rate, 


Y 


whe  re 


k 


r 


-i 


M 

M+2 


(30a) 


(30b) 


and  D  -  (.cT/eBjv .  /SJ  is  the  ion  diffusion  coefficient  based  on  the  total 
X  in  1 

temperature  T  =  Tg  +  The  second  damps  at  the  electron  parallel 

diffusion  rate 


k:- . 

ft 


► 


Y  »  -  k^  D  (31) 

z  le 

V.  GLOBAL  EIGENFUNCTIONS 

In  Section  III  we  derived  a  local  dispersion  equation  [(12)  and  (19)] 
which  was  solved  to  obtain  the  growth  rate  Y(9)  in  Section  IV.  This  local 
growth  rate  is  peaked  on  the  backside  of  the  cloud  (9-0).  We  now 
investigate  under  what  conditions  the  gradient  drift  instability  forms  a 
global  eigenmode  as  it  grows  on  the  cylindrical  waterbag. 


Si 

f; 


22 


j. 


A. 


Such  global  eigenmodea  can  exist  if  solutions  JCs)  can  be  constructed 
which  are  localized  around  some  angle  9q.  To  obtain  these  solutions  we 
take  the  growth  rate  in  (12)  and  (19)  to  be  independent  of  9  so  that  the 
dispersion  equation  yields  k  (9,Y).  We  then  make  the  identification 
k  «  -i  r  ^  3/39  to  obtain  a  differental  equation  for  *(9). 

a  C 

2  2  2 

In  the  limit  X  »  a  k  r  «  1,  (12)  yields  k_.  =•  k^  and  the  local 

2  C  "T  0 

dispersion  equation  in  (19)  becomes 

Yq  3^/39  -  iY  exp(i9)5  +  -  0.  (32) 

Since  this  equation  is  first  order  in  3/39,  there  are  no  bounded  solutions 
so  that  there  are  no  exponentially  growing  global  eigenmodes  when  k^  is 
small.  This  result  is  consistent  with  previous  calculations  in  the  limit 
k^  =•  0  where  it  was  shown  that  the  energy  always  cascades  to  lower  poloidal 
mode  numbers  (Overman  et  al.,  1983).  The  local  growth  rate  [see  (11)]  is 
the  approximate  rate  of  increase  of  the  amplitude  of  a  broad  spectrum  of 
modes  centered  around  k  . 

We  now  consider  the  case  where  X^  >>  1  while  k_  is  sufficiently  small 

so  that  the  terms  proportional  to  v,  /SJ,  in  (12a)  and  (12b)  can  be 

in  1 

2 

neglected.  We  also  assume  that  Y  >>  k  D  and  v  >>  v  , .  With  these 

z  ie  en  e  i 

constraints  (12a)  and  (12b)  simplify  to 


and  (19)  yields  the  equation 

.2  _  a  .  . 

Y^  COS9  ^  (|)  +  [u)„  +  Y  X  3in9]  ^  T  X^  -  0. 


(33) 


We  simplify  this  equation  by  defining  a  new  dependent  variable 


5  »  ()»  exp[5^f(8)] 

with  3f/3e  •  [ojq  *  Yq  ^  3in0]/2YQ  co30  and  find  that 

3^(Ji/30^  +  V(0)({>  -  0  (S'^a) 

where  the  potential  V(0)  is  given  by 

Y  1  “o 

V(e)  -  {(^  0030  ^  5-)  X  -  3in0 

0 

-  [a.„  Y.  X  3in0]^}  — ~  .  (3^b) 

Uy'^  ^  cos  0 

o 

We  first  consider  the  limit  T  -  0  and  expand  V(3)  around  0  »  0, 

V  -  x(^  -  '  F  (35) 

The  bound  state  solutions  for  this  potential  have  eigenvalues 

Y  -  YQ(n+1 ),  (36) 

where  n  is  a  non^negative  integer.  The  mode  is  localized  on  the  backside 

-1  /2 

of  the  cloud  and  has  an  angular  width  d0  -  X  .  <<  1,  which  decreases  with 

increasing  k_.  Thus,  the  expansion  of  V  around  8  -  0  is  valid  for  X  >>  1 . 
Note  also  that  higher  order  modes  n,  which  have  more  structure  in  the 


poloidal  direction,  have  larger  growth  rates.  This  result  is  consistent 
with  local  theory  where  Y  increases  with  the  poloidal  mode  number  m. 

The  potential  V  in  (3^)  is  modified  by  finite  thermal  effects  when 
ujq  -  YqX.  Comparing  the  magnitude  of  the  various  terms  in  (3**)f  we  find 
that  the  term  proportional  to  ^  is  of  order  A  >>  1  larger  than  the 

remaining  terms  unless  the  terms  within  the  bracket  cancel.  We  therefore 
look  for  a  mode  localized  around  the  angle  9q  defined  by 

3in0g  -  -Uq/YqX.  (37) 


Near  0q  the  potential  assumes  the  form 


V(0) 


X 


Y_  1 
Yq  CO30Q 


F  x2(0 


(38) 


The  bounded  solutions  have  eigenvalues 

2 

Y  -  Y  (n>l)  [l  -  (39) 

’o* 

where  n  is  again  a  non-negative  integer.  When  Uq  <<  Y^X,  the  growth 
reduces  to  the  previous  zero  temperature  result  in  (36).  As  Uq  increases 
the  growth  rate  decreases  until  the  mode  becomes  stable  at  uJq  ■  Y^X  or 


1/2 

a 


k  Si 
z  eB 


> 


2 

M+2 


(40) 


The  reduction  of  the  growth  rate  in  (39)  is  a  consequence  of  the 
localization  of  the  mode  at  an  angle  0q  •  0.  The  growth  rate  in  (39)  can 


be  rewritten  as 


Y  -  YqCh+I)  OOS0Q. 

At  the  marginal  stability  point,  8q  -  ir/2,  i.e.,  the  mode  localizes  in  a 
region  where  there  is  no  driving  force.  Above  the  threshold  in  (40)  there 
are  no  bounded  solutions  to  (34)  since  9^  moves  into  the  complex  plane. 

The  physics  behind  the  localization  of  the  mode  can  be  readily 
understood.  As  we  discussed  previously  in  Sec.  IV,  the  convection  of  the 
background  plasma  past  the  cloud  causes  the  mode  to  propagate  with  a 
frequency  [see  (21b)] 


(D  ■  mY-  sine.  ( 41 ) 

r  cj 

On  the  other  hand,  the  thermal  effects  cause  the  mode  to  propagate  with  a 
frequency  [see  (27a)] 


•  iiujIq/X.  (42) 

These  two  frequencies  balance  each  other  at  the  angle  9^  given  in  (37).  The 
mode  localized  at  this  angle  has  no  real  frequency  and  grows  at  the  local 
growth  rate  corresponding  to  this  angle.  Wien  the  propagation  rate  due  to 
the  thermal  effects  is  everywhere  larger  than  the  convection  flow  no 
localized  solutions  exist.  In  this  regime  the  E  »  B  gradient  drift 
instability  is  effectively  stable  since,  unlike  the  limit  where  k^  -  0,  the 
rate  of  propagation  of  the  mode  from  the  unstable  backside  of  the  cloud  to 
the  stable  frontside  exceeds  the  rate  of  growth  of  the  mode.  Convective 
amplification  of  perturbations  in  this  finite  k^  limit  is  therefore  not 


significant. 


We  have  investigated  the  influence  of  finite  parallel  wavelength  on 


the  stability  of  a  cylindrical  plasma  cloud.  We  first  derived  a  dispersion 
relation  for  the  local  growth  rate  y(0)  of  the  E  x  B  gradient  drift 
instability,  where  0  is  the  poloidal  angle  (0  »  C  on  the  "backside”  of  the 
cloud  as  shown  in  Fig.  1).  For  sufficiently  large  values  of  the  parallel 
wavenumber  the  local  growth  rate  is  negative  for  all  values  of  0.  In 
the  cold  plasma  limit. 


(43) 


the  stability  criterion  is  given  in  (25a)  while  in  the  warm  plasma  limit 
Cthe  inequality  in  (43)  is  reversed],  the  stability  condition  is  given  in 
(28b).  In  both  cases  the  lowest  poloidal  mode  numbers  are  stabilized  at 
the  smallest  values  of  k^. 

We  have  also  investigated  under  what  conditions  exponentially  growing 
global  eigenmodes  can  exist.  In  the  limit 

-  ak^r^  «  1  (44) 

z  c 


with  a  =  ^e^i'^'^e'^in  exponentially  growing  solutions.  In  this 

limit  the  energy  cascades  to  lower  poloidal  mode  numbers  as  the  instability 

2 

grows  as  found  by  Overman  et  al.  (1983).  Wien  \  >>  1  a  localized  mode  of 

-1  /2 

angular  width  0  -  X  <<  1  forms  on  the  "backside"  of  the  plasma  cloud 
and  grows  exponentially  in  time.  Thus,  the  cascade  of  energy  to  lower  mode 
numbers  no  longer  takes  place  when  k^  is  sufficiently  large.  The  structure 
of  the  mode  for  this  case  is  illustrated  in  Fig.  6a.  This  figure  is  drawn 


detailed  description 


in  the  rest  frame  of  the  circular  cloud.  The  neutral  wind  is  moving  from 
the  left  to  right  with  a  uniform  velocity  2V^/{M*2)  while  the  background 
plasma  flows  to  the  left  with  a  velocity  The  dashed  line  illustrates 
the  amplitude  of  the  lowest  order  mode.  At  still  larger  values  of  k^  the 
unstable  mode  becomes  localized  at  a  finite  angle  0q  given  by 


where 


u  cT  1/2, 

"3“  “2 


is  the  effective  diamagnetic  velocity  and 


V 


b 


2JL  V 
M+2  n 


with  M  -  n^/njj.  For  a  waterbag  distribution  the  density  scale  length  which 

usually  appears  in  definition  of  the  diamagnetic  drift  velocity  is  replaced 

^1/2 

by  the  radial  scale  length,  a  .  For  or 


1/2.  cT 
a  k  — 
z  eB 


2 

M+2 


(46) 


the  mode  is  completely  stable.  The  physical  mechanism  which  causes  the 
localization  at  0^  as  well  as  the  stabilization  is  illustrated  in  Fig. 
6b.  In  the  local  dispersion  relation  the  diamagnetic  effects  cause  the 
mode  to  propagate  in  the  pololdal  direction  with  a  frequency 


u)  -  k„V  V(M+2) 


where  k.  is  the  poloidal  wavenumber.  The  velocity  of  the  background  plasma 
Just  outside  the  cloud  (in  the  reference  frame  of  the  cloud)  is  given  by 

V  -  V.  sine  e. 

~  b 

The  convection  of  the  background  plasma  past  the  cloud  causes  the  mode  to 
propagate  with  a  frequency 

u)  -  kgVj^  3in0. 

These  two  frequencies  balance  to  produce  a  non-propagating  mode  at  the 
angle  0q  defined  in  (M5).  The  dashed  line  in  Fig.  6b  illustrates  the 
localization  of  the  mode  in  this  case.  When  the  diamagnetic  velocity 
everywhere  exceeds  the  flow  of  the  background  plasma  around  the  cloud  the 
gradient  drift  instability  is  convectively  stabilized. 

Finally,  we  apply  these  results  to  the  structuring  of  barium  clouds 
and  discuss  their  application  to  the  sor^called  "striation  freezing" 
phenomenon  (Linson  and  Meltz,  1972).  Basically,  it  has  been  observed  that 
barium  clouds  released  in  the  ionosphere  structure  because  of  the  E  *  B 
gradient  drift  Instability  and  develop  field-aligned  strlations.  The  first 
generation  of  strlations  can  also  undergo  further  structuring,  at  times, 
and  break  up  into  even  smaller  strlations  (i.e.,  smaller  in  size  transverse 
to  Bq).  This  process,  known  as  bifurcation,  appears  to  continue  until  a 
minimum  transverse  scale  size  is  reached  (which  we  refer  to  as  the 
"freezing  scale  length").  For  barium  clouds  released  at  altitudes  -  180  km 
the  freezing  scale  length  is  roughly  *100  m  (Prettie,  1985).  A  number  of 


studies  have  been  carried  out  which  address  this  problem  (Francis  and 


Perkins,  1975;  McDonald  et  al.,  1981;  Zalesak  et  al.,  1984;  Drake  et  al., 
1985;  Sperling  and  Classman,  1985).  Rather  than  describe  the  detailed 
processes  proposed  in  these  papers,  it  is  sufficient  to  note  that  there  is 
no  generally  accepted  model  of  striation  freezing  at  this  time;  each  model 
has  its  merits  and  shortcomings.  Although  the  theory  presented  in  this 
paper  does  not  explicitly  predict  a  "freezing  scale  length",  we  can  make  an 
estimate  of  this  size  based  on  (46)  and  a  simple  physical  argument.  The 


free  parameter  in  (46)  is  k_ ;  all  other  parameters  are  determined  by 

Z 


ionospheric  conditions.  Thus,  we  need  to  make  a  reasonable  estimate  of 
k^.  We  do  this  by  noting  that  transverse  perturbations  can  map  parallel  to 
the  magnetic  field.  The  relationship  between  parallel  and  perpendicular 


1  /2  1  /2 

scale  lengths  is  approximately  given  by  L  -  (o  /a  )  L  ■  a  L  where  L 

I  I  1  ±  J. 


and  0  refer  to  the  scale  size  and  conductivity,  respectively  (Farley,  1959; 


Goldman  et  al.,  1976).  Assuming  that  k_  -  L  ^  j”"'  we  find  that 


1/2,  ^-1 


(46)  can  be  written  as 


L  < 
± 


cT 

2  eB  V 

n 


(47) 


Thus,  (47)  suggests  that  barium  cloud  strlations  with  transverse  dimensions 

smaller  than  L  would  be  stable  to  further  structuring  by  the  E  *  B  gradn 
L  "  ~ 

lent  drift  Instability.  For  typical  barium  cloud  ionospheric  parameters  at 

180  km,  i.e.,  T  -  Tg  +  T^  -  0.2  eV,  B  -  0.5  G,  -  50  m/sec,  and  M  -  2-10, 

we  find  that  L  -  I60n480  m  which  is  consistent  with  observations.  Of 
X 

1/2  *  1 

course  this  result  is  predicated  on  the  assumption  that  k^  -  [a 

which,  although  plausible,  is  somewhat  ad  hoc.  In  order  to  remove  this 


assumption  it  is  necessary  to  consider  the  finite  length  of  a  barium  cloud 


: VKWB  W:Uni  i.'.  > ' 


which  introduces  a  physical  parameter  which  will  remove  the  arbitrariness 
associated  with  k^.  We  are  presently  developing  such  a  model. 
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DEPT 

01CY  ATTN  DOC  CON  FOR  L-389  R.  OTT 

01CY  ATTN  DOC  CON  FOR  L-31  R.  HAGEI 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS.  NM  87545 

01 CY  ATTN  DOC  CON  FOR  J.  WOLCOTT 

01CY  ATTN  DOC  CON  FOR  R.F.  TASCHEK 

01 CY  ATTN  DOC  CON  FOR  E.  JONES 

01CY  ATTN  DOC  CON  FOR  J.  MALIK 

01CY  ATTN  DOC  CON  FOR  R.  JEFFRIES 

01CY  ATTN  DOC  CON  FOR  J.  ZINN 

01CY  ATTN  DOC  CON  FOR  D.  WESTERVEL' 

01 CY  ATTN  D.  SAPPENFIELO 

LOS  ALAMOS  NATIONAL  LABORATORY 
MS  D438 

LOS  ALAMOS.  NM  87545 
01CY  ATTN  S.P.  GARY 
01CY  ATTN  J.  BOROVSKY 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE.  NM  87115 

01 CY  ATTN  DOC  CON  FOR  W.  BROWN 
01 CY  ATTN  DOC  CON  FOR  A. 
THORNBROUGH 

01CY  ATTN  DOC  CON  FOR  T.  WRIGHT 

01CY  ATTN  DOC  CON  FOR  D.  DAHLGREN 

01CY  ATTN  DOC  CON  FOR  3141 
01 CY  ATTN  DOC  CON  FOR  SPACE  PROJEC 
DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

01CY  ATTN  DOC  CON  FOR  B.  MURPHEY 

01CY  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  DC  20545 

01CY  ATTN  DOC  CON  DR.  YO  SONG 


OTHEH  GOVERNMEHT 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  A  INFO 
AOMIN 

BOULDER,  CO  80303 

Old  ATTN  D.  CROMBIE 
Old  ATTN  L.  BERRY 

NATIONAL  OCEANIC  A  ATMOSPHERIC  AOMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
01 CY  ATTN  R.  GRUBB 
01CY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
Old  ATTN  I.  GARFUNKEL 
01CY  ATTN  T.  SALMI 
Old  ATTN  V.  JOSEPHSON 
Old  ATTN  S.  BOWER 
Old  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

Old  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 

Old  ATTN  L.  SLOAN 
01CY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  9A701 

Old  ATTN  J.  WORKMAN 
Old  ATTN  C.  PRETTIE 
Old  ATTN  S.  BRECHT 

BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 

Old  ATTN  0.  KEISTER 
Old  ATTN  D.  MURRAY 
Old  ATTN  G.  HALL 
01CY  ATTN  J.  KENNEY 


CHARLES  STARK  DRAPER  LABORATORY,  INC. 

555  TECHNOLOGY  SQUARE 

CAMBRIDGE,  MA  02139 
01CY  ATTN  D.B.  COX 
01CY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 

LINTHICUM  ROAD 

CLARKSBURG,  MD  20734 
01CY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OP  ELECTRICAL  ENGINEERING 

ITHACA,  NY  14850 

01CY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
01 CY  ATTN  H.  LOGSTON 
01CY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Wllshire  Blvd. 

Santa  Monica,  CA  90401 
01CY  ATTN  C.B.  GABBARD 
01CY  ATTN  R.  LELEVIER 


ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
01CY  ATTN  J.  ROBERTS 
01CY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORCE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
Old  ATTN  M.H.  BORTNER 
SPACE  SCI  LAB 

GENERAL  ELECTRIC  TECH  SERVICES 
CO.,  INC. 

HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

01CY  ATTN  G.  MILLMAN 


9EOPHXSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01 CY  ATTN  TECHNICAL  LIBRARY 
01CY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVANIA.  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01 CY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN.  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 


KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS.  CO  80933 
01CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET  (P.O  DRAWER  QQ) 
SANTA  BARBARA,  CA  93102 
01CY  ATTN  DASIAC 
01 CY  ATTN  WARREN  S.  KNAPP 
01CY  ATTN  WILLIAM  MCNAMARA 
01 CY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO.  CA  92121 

01CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  A  SPACE  CO..  INC 
P.O.  BOX  504 
SUNNYVALE.  CA  94088 
01CY  ATTN  DEPT  60-12 
01CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  A  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO.  CA  94304 


ALEXANDRIA. 

VA  22311 

OICY 

ATTN 

MARTIN  WALT  DEPT  52-12 

OICY 

ATTN 

J.N.  AEIN 

OICY 

ATTN 

W.L.  IMHOF  DEPT  52-12 

h! 

01CY 

ATTN 

ERNEST  BAUER 

OICY 

ATTN 

RICHARD  G.  JOHNSON 

OICY 

ATTN 

HANS  WOLFARD 

DEPT 

52-12 

■ 

OICY 

ATTN 

JOEL  BENGSTON 

OICY 

ATTN 

J.B.  CLADIS  DEPT  52-12 

INTL  TEL  A  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

01CY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01 CY  ATTN  THOMAS  POTEMRA 
01 CY  ATTN  JOHN  DASSOULAS 


MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO.  FL  32805 

01 CY  ATTN  R.  HEFFNER 

MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CY  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOULE 
01CY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
01CY  ATTN  R.W.  HALPRIN 
01CY  ATTN  TECHNICAL 

LIBRARY  SERVICES 
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IflSSIOM  RISEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA.  CA  93101 
01CI  ATTN  P.  FISCHER 
OICT  ATTN  H.P.  CREVIER 
01 CT  ATTN  STEVEN  L.  CUTSCHE 
01 CY  ATTN  R.  BOCUSCH 
OICY  ATTN  R.  HENDRICK 
01 CY  ATTN  RALPH  KILB 
OICY  ATTN  DAVE  SOWLE 
OICY  ATTN  F.  FAJEN 
OICY  ATTN  M.  SCHEIBE 
OICY  ATTN  CONRAD  L.  LONGHIRE 
OICY  ATTN  B.  WHITE 
OICY  ATTN  R.  STAGAT 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD.  S.E. 
ALBUQUERQUE.  NM  87106 
OICY  R.  STELLINGHERF 
OICY  N.  ALME 
OICY  L.  WRIGHT 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN.  VA  22101 
OICY  ATTN  W.  HALL 
OICY  ATTN  W;  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES.  CA  9002S 

OICY  ATTN  E.C.  FIELD.  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK.  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH 

PHOTOMETRICS.  INC. 

4  ARROW  DRIVE 
WOBURN.  MA  01801 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS.  INC. 

P.O.  BOX  3027 
BELLEVUE.  WA  98009 

OICY  ATTN  E.J.  FREMOUW 


PHYSICAL  DYNAMICS.  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

RAD  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REX,  CA  90291 
OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  WILLIAM  B.  WRIGHT,  JR 
OICY  ATTN  WILLIAM  J.  KARZAS 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 

RAND  CORPORATION.  THE 
15450  COHASSET  STREET 
VAN  NUYS,  CA  91406 

OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY.  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK.  NY  10036 

OICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA.  CA  92037 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  E.  FRIEMAN 
OICY  ATTN  E.A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS,  INC 
LAB  1710  GOODRIDGE  DR. 

MCLEAN.  VA  22102 
OICY  J.  COCKAYNE 
OICY  E.  HYMAN 


SJtl  INTERMATIONAL 
333  RAVENSWOOD  AVENUE 
HENLO  PARK.  CA  9A025 
01CY  ATTN  J.  CASPER 
01 CT  ATTN  DONALD  NEILSON 
01CT  ATTN  ALAN  BURNS 
01 CX  ATTN  G.  SHITH 
01 CX  ATTN  R.  TSUNODA 
01 CX  ATTN  DAVID  A.  JOHNSON 
01 CX  ATTN  HALTER  G.  CHESNUT 
01CX  ATTN  CHARLES  L.  RZNO 
01 CX  ATTN  HALTER  JAXE 
01CX  ATTN  J.  VZCKREX 
01CX  ATTN  RAX  L.  LEADABRAND 
01CX  ATTN  G.  CARPENTER 
01CX  ATTN  G.  PRICE 
01CX  ATTN  R.  LIVINGSTON 
01CX  ATTN  V.  GONZALES 
01 CX  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  HIGGINS  AVENUE 
BEDFORD.  MA  01730 

01CX  ATTN  H.P.  BOQUIST 

TRH  DEFENSE  A  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH.  CA  90278 
01CX  ATTN  R.  K.  PLEBUCH 
01CX  ATTN  S.  ALTSCHULER 
01 CX  ATTN  0.  DEE 
01 CX  ATTN  D/  STOCKHELL 
SNTF/1575 

VISIDXNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MA  01803 
01CX  ATTN  H.  REIDX 
01CX  ATTN  J.  CARPENTER 
01 CX  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

01 CX  ATTN:  N.  ZABUSKX 

DTIC 

02CY 

CODE  1220 
OICY 

DIRECTOR  OF  RESEARCH 
U.S.  NAVAL  ACADEMY 
ANNAPOUS,MD  21402 
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